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We performed high-dimensional quantum dynamical cal- 
culations of the dissociative adsorption and associative des- 
orption of H2/Cu(lll). The potential energy surface (PES) 
is obtained from density functional theory calculations. Two 
regimes of dynamics are found, at low energies sticking is de- 
termined by the minimum energy barrier, at high energies by 
the distribution of barrier heights. Experimental results are 
well-reproduced qualitatively, but some quantitative discrep- 
ancies are identified as well. 
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Dissociative adsorption is one of the crucial steps in 
heterogenous catalysis. This has motivated, apart from 
its fundamental interest for surface chemistry, a large 
number of studies on well defined model systems where 
the dissociative adsorption process can be investigated in 
detail. The system which has been studied to the largest 
extent is H2/Cu(lll) (see, e.g., ||,^). For this system 
the adsorption is hindered by a noticable energy bar- 
rier. It has been found that the dissociative adsorption 
of hydrogen approximately obeys normal-energy scaling 
Q , which means that the sticking probability is only de- 
termined by the component of the molecule's momen- 
tum perpendicular to the surface. Usually normal-energy 
scaling is associated with a flat, structureless surface. For 
the (111) metal surface this seemed to be an appropriate 
description as this is the most densely packed surface. 

Indeed the electron density at clean metal (111) sur- 
faces shows only a small corrugation of less than 0.1 A at 
a distance of 1 A above the center of the surface layer [Q . 
However, adsorption properties are determined by the 
local chemical interaction on the surface. In recent to- 
tal energy calculations [pj using density functional theory 
(DFT) together with the generalized gradient approxima- 
tion (GGA) for the exchange-correlation functional, 
the potential energy surface (PES) for H2/Cu(lll) has 
been determined. Besides the GGA the main approx- 
imation is due to the super-cell approach and the k- 
summation. A rather strong corrugation is found when 
all six degrees of freedom of the II2 molecule are var- 
ied. For a 4-layer substrate with the molecular axis kept 
parallel to the surface, which is the most favorable config- 
uration for dissociation, the energy barrier to adsorption 
varies between 0.73 eV and 1.43 eV within the surface 
unit cell. Also the distance of the barrier from the sur- 



face and the extension of the hydrogen bond at the bar- 
rier region strongly depend on the point of impact in the 
unit cell and the orientation of the molecule. 

In order to investigate the consequences of the strong 
corrugation on the adsorption and desorption dynamics, 
we have taken the GGA-PES as an input for a quan- 
tum dynamical simulation. We use a five-dimensional 
parametrization of the results obtained for the 4-layer 
substrate with the molecular axis being kept parallel to 
the surface as for this configuration the lateral corru- 
gation of the PES has been determined in great detail 
in the DFT calculations. Although the polar rotation 
of the molecule is not considered, to our knowledge this 
is the first study of the influence of lateral corrugation 
on the quantum dynamics of adsorption and desorption 
which takes into account such a large number of degrees 
of freedom using a realistic PES. Because of the high- 
dimensionality of the calculations, the comparison of the 
results with experimental data for the Il2/Cu(lll) sys- 
tem will in addition provide information about the accu- 
racy of the theoretical GGA-PES. We will show that the 
GGA-PES reproduces the experimental findings qualita- 
tively well, but some quantitative differences are identi- 
fied. We will demonstrate the importance of the corruga- 
tion of the PES for the adsorption dynamics, the angular 
dependence of the sticking probability, the mean kinetic 
energy in desorption, and the desorption flux. 

The quantum dynamics is determined by solving 
the time-independent Schrodinger equation for the II2 
molecule moving in the GGA-PES. We use the concept 
of the local reflection matrix (LORE) H and the inverse 
local transmission matrix (INTRA) |^. The propaga- 
tion of the reflection and inverse transmission matrices 
on an accurate piecewise constant representation of the 
PES avoids exponentially increasing evanescent waves 
whereby numerical instabilities are prevented. The cen- 
ter of mass distance of the molecule from the surface Z, 
its two surface coordinates X and Y and the interatomic 
molecular spacing r are taken as dynamical variables. 
The molecular azimuthal orientation is treated as a fixed 
parameter in the adsorption simulation, and all theoret- 
ical results are then obtained by averaging over twelve 
different orientations. This is a reasonable approxima- 
tion because the molecular rotation is a slow process. 
For instance, an II2 molecules with a kinetic energy of 
0.5 eV in the rotational state J = 4 corresponding to 
an rotational energy of 0.15 eV rotates only once in the 



time it takes to move forward 4.2 A, while the interaction 
range of the potential is 0.5 A. 

The local nature of LORE and INTRA does not only 
give information about final transition probabilities, but 
also about the scattering solution of the wave function 
everywhere in space 0. To give an idea about the size 
of the corrugation of our PES and its consequences on 
the dynamics we have plotted in Fig. [^ the positive real 
part of the wave function for an incident beam with total 
kinetic energy Ei = 1.08 eV and incident angle 9i — 45° 
scattered at a surface with one-dimensional lateral corru- 
gation having the same variation in energy and location 
of the barrier as the full potential. The sticking proba- 
bility is given by the flux of the particles that traverse 
the barrier. To make both the incoming plus transmitted 
and the reflected part of the beam visible in one single 
plot, we have restricted in the plot the lateral extension 
of the beam perpendicular to its propagation direction, 
which is in principle infinite. The splitting of the beam 
in a refiected and a transmitted part can be traced very 
clearly. Also the development of the different diffraction 
beams can be followed. 

In Fig. we show results for the sticking probability 
as a function of the normal kinetic energy of the inci- 
dent II2 beam initially in the vibrational ground state 
for various polar angles of incidence {9i = 0° refers 
to normal incidence). Furthermore we present results 
for two-dimensional calculations corresponding to non- 
corrugated surfaces with the minimum and maximum 
barrier of the fully corrugated surface. Also a stick- 
ing curve derived from experiment llfl] is shown. For 
9i = 0° two energy regimes can be separated. At low 
energies (Fig. 0b) the "5D" sticking curve for the corru- 
gated surface is parallel to the 2D sticking curve for the 
non-corrugated surface with the minimum energy barrier. 
The 2D sticking curve for the maximum barrier is sup- 
pressed by more than ten orders of magnitude and has a 
different slope. These results indicate that the "5D" data 
can be explained by a kind of key-hole effect. Since stick- 
ing at these energies can only be achieved by tunneling, 
it is very sensitive to the barrier height. The sticking for 
normal incidence at low energies is therefore determined 
by the molecules which hit the surface within the low- 
est energy configuration; molecules which hit the surface 
at unfavorable sites or with unfavorable orientation do 
practically not contribute to the sticking. 

For incident energies larger than approximately 0.6 eV 
the slope of the sticking curve for the corrugated sur- 
face in the logarithmic plot changes. At this energy ad- 
sorption without tunneling becomes possible since the 
effective minimum barrier height is reduced by 0.15 eV 
for molecules in their vibrational ground state due to 
the softening of the II2 bond at the surface and the ac- 
companied reduction of the zero-point vibrational energy. 
In the linear plot (Fig. 0a) the sticking curve obtained 
for the corrugated surface connects the onset of stick- 



ing for the minimum barrier approximately linearly with 
the point where the sticking probability in the maximum 
barrier system reaches unity. Test calculations with a 
changed barrier widths but unchanged barrier heights 
yield almost unchanged sticking probabilities at these en- 
ergies. These facts indicate that the sticking probability 
in this energy range is determined by the distribution of 
the barrier heights. This is the first multidimensional 
evidence that for large energies the sticking can be un- 
derstood in terms of the region of the surface that clas- 
sically is available to dissociation, which is an important 
assumption in the so-called hole-model Q . 

The comparison of the experimentally derived stick- 
ing curve with the theoretical results for the high- 
dimensional PES shows that the increase for sticking 
probabilities between 0.01 and 0.4 is actually similiar. 
The main difference in this range of sticking probabilities 
is that the theoretical results appear at energies which are 
too high by about 0.2 eV. This indicates that the min- 
imum energy barrier for dissociation of the GGA-PES 
should be lowered from 0.73 eV to about 0.5 eV. In fact, 
convergence tests reported in Ref. M have shown that 
such a lowering of the barrier is to be expected. 

For low and for high energies, however, the theoret- 
ical sticking probabilities are at variance with the ex- 
periment. The experimental curve has been derived by 
fitting existing adsorption data and desorption data via 
the principle of detailed balance to an assumed func- 
tional form of the sticking curve |jlQ]. In detailed bal- 
ance the desorption flux is determined by the product of 
the sticking probability, which increases exponentially in 
the low-energy regime, with the exponentially decreasing 
Boltzmann factor, resulting in a desorption flux which is 
strongly peaked at energies close to the minimum barrier 
height. Beam adsorption data, which are only available 
up to energies of about 0.5 eV for H2 |10|, are domi- 
nated at low energies by the initially vibrationally excited 
molecules [g[. Thus at least for initially non- vibrating 
molecules the fitting procedure might not be very sen- 
sitive at low and high energies. For energies close to 
the minimum barrier height the agreement between our 
results and experiment (except for the energy shift) sug- 
gests that the spatial corrugation of our PES resembles 
that of reality. This also implies that the width of the ex- 
perimentally derived sticking curve does not correspond 
to the barrier width ||l^ , but to the distribution of bar- 
rier heights since it is this distribution which determines 
the increase of the sticking probability in this range. 

Just recently it has been pointed out |l3] that, apart 
from the change in momentum of inertia, an isomorphism 
exists between surface corrugation and a planar rotor. 
This can be used to give an estimate about the effect of 
dynamically including the rotations. The total energy 
for a rotated H2 dimer has been calculated: At the opti- 
mum reaction configuration the energy required to turn 
the molecule upright is 0.4 eV S which is of the or- 



der of the lateral corrugation. Thus taking into account 
the rotation will have a similiar effect as including the 
surface corrugation, namely that for J = in the low 
energy regime the sticking probability will basically be 
suppressed by a constant factor due to the smaller key- 
hole in phase space, while in the high energy regime the 
energy at which the sticking probability reaches one will 
be shifted to larger energies due to the existence of higher 
barriers. The onset of the high energy regime will not be 
changed since the minimum barrier will remain the same. 

For non-normal incidence also two energy regimes can 
be separated. At high energies (above 0.6 eV, see Fig. ^a) 
the sticking obeys normal-energy scaling in agreement 
with the experiment although the surface is strongly cor- 
rugated. This issue has already been addressed by Dar- 
ling and HoUoway |lj]. Using model potentials they 
showed that the variation of the barrier height and the 
variation of the barrier position have opposing effects as 
far as the role of addional parallel momentum for ad- 
sorption is concerned. If the maximum barriers are at a 
greater distance from the surface these opposing effects 
cancel to a large extent leading to approximate normal- 
energy scaling even for strongly corrugated surfaces [Q . 
The GGA-PES shows these features (see Fig. |l|) indicat- 
ing that this explanation also applies to our calculations 
including the full lateral corrugation. 

However, Fig. 0b reveals that in the low energy range, 
which had not been considered in Ref. |14|, additional 
parallel momentum helps to cross the barrier. The situ- 
ation is shown in Fig. ^. A large part of the transmitted 
fraction of the beam propagates through the minimum 
barrier site without moving like a classical particle, but 
being steered in the tunneling process to the energeti- 
cally most favorable propagation path thereby convert- 
ing parallel into normal energy, which in turn increases 
the tunneling probability ||l5|l . 

In Fig. y the mean kinetic energy of desorbing 
molecules versus desorption angle at a surface temper- 
ature of Ts = 925 K is plotted. Measurements by Comsa 
and David |16[| have shown that the kinetic energy is inde- 
pendent of the desorption angle, which is surprising tak- 
ing into account normal energy scaling for Il2/Cu(lll) 
in a classical model (solid line in Fig. g). This issue 
has been discussed at length by Michelsen and Auer- 
bach |1^] who showed via detailed-balance arguments 
that due to the significant width of the sticking prob- 
ability, which effectively changes with desorption angle, 
a constant mean kinetic energy in desorption can be 
achieved even under the assumption of normal-energy 
scaling. Our calculations also give a kinetic energy of 
desorbing molecules almost independent of the desorp- 
tion angle. Here it is mainly due to the inclusion of 
geometric corrugation which makes parallel momentum 
more effective for traversing the barrier leading to a more 
isotropic distribution of desorbing molecules |15[ . 

The mean kinetic energy obtained in our "5D" quan- 



tum dynamical calculation is 0.2 eV larger than the one 
measured in experiments |lq,0| . This again corresponds 
to the shift in the minimum energy barrier height of our 
PES as discussed above. The vibrational excitation of 
desorbing molecules at T^ = 925 K is calculated to be 
0.87% while the experimental result is 2.9% 0. The 
GGA places the energy barrier for dissociation correctly 
in the "exit" channel in the rZ-plane (see Fig. la in M). 
However, the H2 bond length at the barrier top is only 
extended by 33% so that the barrier position is still too 
early |l|jl|]. 

Calculations using only a two-dimensional PES with 
the minimum energy path yield 0.67 eV for the mean 
kinetic energy and 0.83 % for the vibrational excitation 
in desorption at Ts = 925 K. These numbers are rel- 
atively close to the ones derived for a fully corrugated 
surface indicating that most desorbing molecules follow 
the minimum energy path. This shows that it is justi- 
fied to neglect corrugation in the simulation of desorption 
processes if there is one preferred reaction path to des- 
orption which is energetically much more favorable than 
all others. 

The calculated angular distribution of desorbing 
molecules can be fitted to a cos" 9-cuTve with n « 25. 
This distribution is narrower than the one seen in exper- 



iments which yield an exponent of n^ 



8|iq. This can 



partly be attributed to the fact that the kinetic energy 
of desorbing molecules is too high in our calculations. 
Furthermore, experiments have shown that the angular 
distribution in desorption is broadened if the motion of 
the surface atoms is taken into account [E0|, which is 
neglected in our calculations. 

In conclusion, we have presented a dynamical study 
of the dissociative adsorption and associative desorption 
in the system H2/Cu(lll) using a five-dimensional po- 
tential energy surface which has been determined by ab 
initio calculations in the generalized gradient approxi- 
mation. Except for an energy shift of 0.2 eV due to 
the restricted super cell and k-point set, the GGA-PES 
describes the surface corrugation qualitatively well. In 
the tunneling regime sticking is determined by the min- 
imum energy barrier. For energies larger than the ef- 
fective minimum barrier height the width of the sticking 
curve corresponds to the distribution of barrier heights, 
and sticking obeys normal-energy scaling although the 
PES is strongly corrugated. The mean kinetic energy of 
desorbing molecules is almost independent of the desorp- 
tion angle due to the variation of the barrier location. 
In order to simulate desorption experiments corrugation 
can be neglected if there is one energetically favorable 
reaction path to desorption. 
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FIG. 3. Mean kinetic energy of desorbing molecules versus 
desorption angle at a surface temperature of Ts = 925 K 
obtained in our quantum dynamical calculation (open circles) 
and in a classical model for a flat surface (after Ref. |lq]) with 
the minimum energy barrier of 0.73 eV. 



FIG. 1. Positive real part of the wave function of H2 scat- 
tered at a PES with one-dimensional lateral corrugation sim- 
ulating H2/Cu(lll). The incident beam has a kinetic energy 
of Ei = 1.08 eV, the incident angle is 9i = 45°. The contour 
spacing for the potential (thick lines) is 0.1 eV. The positions 
of the Cu atoms are at the potential maxima. 



FIG. 2. Sticking probability versus normal kinetic energy 
for molecules initially in the vibrational ground state, a) 
Linear plot, b) logarithmic plot (note the different energy 
range). 5D-calculations for different incident angles at the 
corrugated surface: solid line Oi ^ 0°, O Oi = 15° , A Oi = 30°, 
O di = 45°; 2D-calculations corresponding to a flat surface: 
dcish-dotted line for the minimum barrier, long-dashed line for 
the maximum barrier. Dashed line: Experimental curve in- 
corporating adsorption data and, via the principle of detailed 
balance, also desorption data (from Ref. 
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